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0\ : ABSTRACT 

Qh' We derive the equation for the surface of equal arrival time of radiation from a 

. thin relativistic shell interacting with an external medium, representing the afterglow 

i of a gamma-ray burst produced by a fireball. Due to the deceleration, these surfaces 

CN i become distorted ellipsoids and, at sufficiently late times, most of the light (either 
bolometric or in a given band) comes from a ring-like region whose width depends 

J> ■ only on age. We analyze the shape of these surfaces and the radiation received 

■ from different angles for different dynamic and radiative regimes and homogeneous 

■ or power-law external densities. We calculate angle-integrated bolometric and fixed 
^ . frequency fluxes, and we tabulate the most relevant parameters that describe the equal 

, arrival time surfaces and the source brightness distribution, quantities that are useful 

0^ 



for more accurate analytic estimates of the afterglow evolution. 
O I Subject headings: gamma-rays: bursts - methods: analytical 



1. Introduction 



^ I The afterglows of gamma-ray bursts (GRB) appear to be well-fitted by decelerating relativistic 

fireball models (Tavani, 1997; Vietri, 1997; Waxman, 1997a; Wijers, Rees & Meszaros , 1997). 
This picture (Meszaros & Rees, 1997), in its simplest form, assumes that the bulk of the radiation 
comes from the external blast wave pushed ahead of the fireball with a diminishing bulk Lorentz 
factor, which is predicted to produce radiation at wavelengths longer than 7-rays decaying as 
a power-law in time, in good agreement with observations. Two interesting consequences of 
the deceleration dynamics are that most of the late radiation comes from a narrow ring, rather 
than the entire visible surface (Waxman, 1997b), and that the usual estimate for the transverse 
size of a relativistically expanding cloud under-estimates the real one (as we show in §^). This 
has consequences for the apparent expansion rate of the fireball, the evolution of scintillation 
properties of the radio-emitting remnant (Goodman, 1997; Frail et al, 1997), and the probability 
of microlensing of GRB afterglows (Loeb & Perna, 1997). The exact shape of the surface, the 
dimensions and the expansion rate depend on the dynamic and radiative regimes, as well as 
on assumed properties of the external medium. We present detailed analytical and numerical 
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calculations of these properties taking into account the effects of the dynamic and radiative 
efficiency regimes, for both homogeneous and power-law density external media. We investigate 
the equal arrival time surfaces and source apparent width evolution, provide estimates of the 
width of the rings, and present simple analytic expressions for both the "average" longitudinal 
and transverse sizes, in the cases of either bolometric or fixed frequency band observations. 



2. Equal Arrival Time Surfaces 

Here we derive the equation of the surface that the observer sees at given time T. For 
simplicity, we assume that the source of radiation can be approximated as a surface (see §Q for a 
discussion on this approximation). We also assume the external medium to be isotropic, but not 
necessarily homogeneous; therefore at any lab-frame time t (measured in the center of explosion 
frame), the ejecta is spherical. The observer equal-T surface is symmetric with respect to the line 
of sight (l.o.s. ) toward the center of explosion which released the relativistic ejecta, therefore its 
equation is given by two coordinates: a polar angle 9 measured from this central l.o.s. and a radial 
coordinate r. In the absence of deceleration, the equal-T surface is an ellipsoid (Rees 1966) with 
semi-major axis T'^[3cT and semi-minor axis T(5cT, where F = (1 — is the constant Lorentz 

factor of the freely expanding ejecta and c is the speed of light. For large F, the ellipsoid is very 
elongated (high eccentricity). When deceleration is present, the shape of the equal-T surface 
departs from that of an ellipsoid. 

The equation of the equal detector time T surface is 

ct -r{t) cos e = cT , (1) 

where r and t are related by dr = [ighcdt, with /J^/jC the speed of the shock that sweeps up the 
external medium. The Lorentz factor 7s/j of this shock can be approximated (e.g. Meszaros , Rees 
& Wijers, 1997) as a power-law in r: 

lsh = Tsh{r/raec)-'' ; ^=^>0> (2) 

1 + 

where a (< 3) is the index of external gas density power law dependence (p oc r~"), 5 describes 
the dynamics (5 = for momentum-conserving and 5 = 1 for energy-conserving evolution), 
and rdec is the deceleration radius, defined as the radius at which the ejecta has swept up a 
mass equal to a fraction Fq ^ of its own mass, Fq being the initial Lorentz factor of the ejecta. 
Numerical simulations of the hydrodynamic interaction between the ejecta and the external 
medium (Panaitescu & Meszaros , 1997) show that ^sh slowly decreases below the deceleration 
radius rdec and that Tsh — (2/3)Fo. Thus equation (^) is correct only for r > rdec- 

For relativistic shocks, the relationship between t and r has a simple form: 
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Substituting t from equation (^) in equation (|l|), the equal-T surface is 



1^2r,;,V « 2n+lj ' ^ ' 

where the reduced variables a = r/r^ec and r = T/T^jec with T^e^ = rrfec/(2r^^c) have been used. 
Tdec is the time-scale for the onset of the deceleration given in equation At given T, the fluid 
moving directly towards the observer [9 = 0) is located at 

x^a. = [(2n + l)T]i/(2"+i)rrf,e, (5) 

this being where the radius is largest and Lorentz factor is smallest on the surface, 
lsh,o = [(2n + 1) r]""^*-^""''^^ r^/i. The relativistic shock condition 7^/^ > 2 used in deriving 
equation (|3|) limits the applicability of equation (^) to r < (2n + l)~^(rs/i/2)^+^/". As an example, 
we show in Figure 1 the equal-T surfaces at different values of r, for To = 500 (corresponding 
to Tgh = 330), for a homogeneous external medium (a = 0). The times indicated in the legend 
correspond to 3.6 hours, 1.5 days, 5.0 days (and 15 days, right panel), if Tdec is taken as 13 
seconds, which would be the value in a burst arising from Tq = 500, an initial energy release of 
£"0 = 10^^ ergs/sr, and a redshift z = \. (A jet of angle 6 jet reduces the energy required by a factor 
^Jet without changing the results, as long as 7^/1 ^ ^Jet)- 

It is customary in analytic derivations to consider that at a given time T the emitting 
surface is located at r = 27^ (T)cT, and that the disk seen by the observer has a radius 
R = r/{2^{T)) = ^{T)cT, as it would be in the absence of the deceleration (i.e. an ellipsoid), 
and to calculate the properties of the received radiation using the physical parameters (magnetic 
field, electron and flow Lorentz factor etc.) of the fluid located at (x = r, y = 0), the center of the 
projected surface. When deceleration is present, the radial coordinate Xmax of the center of the 
equal-T surface can be related to T by integrating dT = dr/(27^^c) and using equation (^): 

= 4(2n + 1)72 cT , (6) 

where we used the flow Lorentz factor 70 = j{xmax) instead of that of the shock 7^/1,0 = V2jo- 
Therefore Xmax is larger by a factor 2(2n + 1) then the value that is typically used, 2'y'^cT. For 
n = 1.5 (adiabatic remnant and homogeneous external gas) one obtains in equation (^) a factor 
16 (Sari 1997), but this factor could be as large as 28 for a strong coupling radiative remnant 
(n = 3). This differs substantially from the usual factor 2 in the longitudinal size. To estimate the 
corresponding departures in the largest transverse size Umax observable assuming the geometry of 
an ellipsoid for the equal-T surface in a decelerating fireball, we consider the projection of this 
surface on the plane perpendicular to the central l.o.s. 

yZx = Xmax/{2io) = 2{2n + l)7ocT . (7) 

The numerical factor in equation (|^) is 8 for n = 1.5. Waxman (1997b) argued that such a 
large transverse size is incompatible with observations and that the correct transverse source 
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size is smaller by a factor 4 than that inferred using a factor 16 in equation (^). The large 
factor in equation (^) is due to the inappropriate use of the geometry of an ellipsoid in the case 
of a decelerated expansion. Since the transverse source size is important for self-absorption 
considerations, for the evolution of the timescale and amplitude of afterglow radio scintillations 
and for gravitational microlensing, it is worth calculating accurately the coefficient in equation (|^ 
that gives the source size. 



3. Characteristic Radiation and Dimensions 

In order to determine the properties of the radiation received by the observer, one has to 
integrate over the surface of equal arrival time the emission from different parts of the shocked 
fluid, taking into account relativistic effects and the fact that each infinitesimal ring [6, + dO] is 
characterized by different physical parameters (magnetic field, electron density, electron Lorentz 
factor, flow Lorentz factor). We assume that the electrons cool only through synchrotron radiation 
(our numerical simulations show that this is a good approximation), and that they are either 
in the radiative or in the adiabatic regime. In the former case we take into account that the 
remnant can be either radiative or adiabatic, depending on the strength of the coupling between 
the radiating electrons, protons and the magnetic field (for more details on the radiative regime 
and the dynamics, see Meszaros , Rees &; Wijers, 1997), while in the latter case the remnant can 
be only adiabatic (energy conservation). Using the scaling relationships for the magnetic field 
B' oc ^r""/^, the co- moving electron density n'^ oc jr~", synchrotron cooling time-scale t'^y oc rj~^, 
expansion time-scale t'^^p oc 7~^r, synchrotron power P^y oc 7^r~" and peak of synchrotron 
spectrum i/'p oc ^-^r^"/^, where 7 is the flow Lorentz factor, the co-moving spectral intensity at the 
synchrotron peak oc n'^{P'sy / I'p) 

min {t^j^, fg^p} is 

^, I 7^ir^"/^ for radiative electrons 
1 ^^i-3q/2 £qj, adiabatic electrons 

The observed spectral intensity at the detector frame peak of the synchrotron spectrum 
= ^pl [7 (1 - /? cos 6)] is ly^ = {up/u'pfll^, therefore 

''^ I a-(4«+i-5-i)(l + 7^0^)-3 t',y > ti^p ■ 

The bolometric co-moving intensity is /' ~ Ipi^'p, and the observed bolometric intensity is 
/ {vp/Vp)^I' . The synchrotron spectrum is approximated as a broken power-law: 1^ = (u/upYl^^ 
below the peak [v < Up) and ly = {v/vp)~'^Iy^ above the peak {v > Up). We considered that 
e = — 1/2 and e = p/2 for radiative electrons and e = 1/3 and e = (p — l)/2 for adiabatic electrons, 
p being the index of the power-law distribution of electrons (we used p = 2.5). The detector flux 
density is obtained by integrating 1^ over the equal arrival time surface, taking into account the 
solid angle subtended by each infinitesimal ring [9, 9 + d9] on this surface. 
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We study first the bolometric brightness distribution on the equal-T surface. Figure 1 shows 
where most of the radiation comes from: the upper half highhghted zone radiates 50% of the total 
energy; 25% of it is emitted by the cap extending from = up to the indicated region and the 
other 25% is radiated by the area extending toward the origin. Similarly, the lower half highlighted 
part radiates 80% of the energy received at detector. Projecting any of these zones on the plane 
perpendicular to the central l.o.s. results in rings that are thin compared to the disk of the entire 
projected surface seen by the observer. 

Apart from the errors mentioned at the end of §§, a new source of inaccuracy in analytic 
calculations of the source size and the radiation received comes from the fact that the observer 
does not receive most of the flux from the central l.o.s. of the deformed ellipsoid. There is a 
significant difference between the average radial coordinates of the regions highlighted in Figure 1 
and that of the fluid on the l.o.s. toward the center, and one must keep in mind that all relevant 
radiation parameters are power-laws in 7, which is a power-law in r. We can assess the error in 
equation (|^), obtained using the "ellipsoid approximation", by calculating an intensity- weighted 
average transverse coordinate y on the equal-T surface and comparing it to jocT. We also calculate 
an intensity-weighted average longitudinal x coordinate on the same surface and compare it with 
Xmax from equation (P), so that the factor (xmax/'S:)^ will estimate the difference between an 
average 7 that should in general be used instead of 70. These averages are given in Table 1, which 
also lists the width w of the outer ring of the source projection on the plane perpendicular to 
the central l.o.s. , which contains 50% of the entire flux received at detector. Also in Table 1 are 
the ratios between the intensity-averaged synchrotron peak frequency 17^ over the T-surface, and 
1^0 = Vp{G = 0). We find that the coefficient 2(2n -|- 1) in equation (^) over-estimates the true value 
of the transverse size y by a factor 2.1 — 2.7 . 

We consider now observations in a fixed frequency band. Figure 2 (upper graphs) shows 
the transverse distribution of the observed synchrotron peak frequency Vp for various constant-T 
surfaces. The unknown coefficient in the expression for Up was chosen so that at T ~ 1 day (with 
Tdec = 13 s), most of the energy the observer receives is in the optical range (~ 2 eV). The lower 
graphs show the peak (or bolometric) luminosity integrated up to the transverse coordinate y 
as a function of y. The ring is indicated by the steep rise in the integrated luminosity, during 
which Vp varies by approximately one order of magnitude around the peak frequency Vpiymax) of 
the radiation from the region that is seen tangentially by the observer. If observations are made 
at energies ^ 10~^fp(ymax) (e.g. in radio, for the times chosen in Figure 2) then the observer 
practically sees only the low-energy part of slope -1/2 (or 1/3) of the synchrotron spectrum 
of the radiation emitted from most points on the equal-T surface, and the entire disk appears 
almost equally bright. However, if observations are made at energies ^ ^Qvpijjmax) (optical or 
X-ray for Figure 2), then the observer sees mainly the high-energy tail of slope —{p— l)/2 (or 
—p/2) of the synchrotron spectrum from the power-law distribution of electrons, and the visible 
region reduces to a ring. For a given observed frequency band, as the shocked fluid is decelerated, 
i^piUmax) crosses the observed band, and the region radiating in that band shrinks from the full 
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disk to a narrow ring with outer boundary at Umax, the edge of the radiating surface. During 
this transition x/xmax decreases while y/umax increases. At energies far above or below i'p{ymax)-, 
these quantities and the width of the "visible" zone are approximately constant in time. Table 2 
gives the asymptotic range of the same coefficients as Table 1, for observations made at a given 
frequency. The first number in each column gives the value of the coefficient when the source 
is seen as disk {u <^V^ ^ t^piVmax), larger width w), and the last number gives the asymptotic 
value of the coefficient when the source has reduced to a ring (v S> I'piymax), smaller w). The 
coefficients have the same range for all frequencies. The particular frequency of the observing band 
only determines the time when the gradual transition between the two coefficients is made, earlier 
in X-rays (few hours) than in optical (~ 1 day) or radio (~ 10 days). It can be seen that the 
radiative remnant gives narrower rings, and that the ring is wider for expansion into a decreasing 
density medium (e.g. a = 2) than into a homogeneous medium (a = 0). 

4. Discussion 

The main conclusions to be drawn from the calculations presented here are: 

1. For the afterglow of a fireball, the equal arrival time surfaces are distorted ellipsoids whose 
shape and evolution depend on the dynamical regime of the remnant, the electron radiative 
efficiency and the density distribution of the external medium. 

2. Equation (|6|) should be used to relate the observer time T with the radial coordinate of the 
center of this surface, and the size of the source should be calculated with an equation similar 
to (0) but with the coefficient 2(2n + 1) replaced by the coefficients in Table 1 for bolometric 
observations. For band observations these coefficients change in time; Table 2 gives only their 
upper and lower limits. This is of relevance for the scintillation of the afterglow in radio (Goodman, 
1997; Frail et al, 1997). 

3. The spectrum and intensity should be calculated using the gas parameters characteristic of the 
ring (rather than those of the l.o.s. to the center). The mean ratio between average peak frequency 
in the ring and in the central l.o.s. is given in Table 1. For a given dynamic and radiative regime 
this ratio is constant in time, thus the power-law time dependences of the observable flux predicted 
by the fireball afterglow models (Meszaros & Rees, 1997; Meszaros , Rees & Wijers, 1997) remain 
unchanged. 

4. When observations are restricted to a narrow energy band, the shape (ring or disk) of the source 
seen by the observer is dependent on the observational band. In any band, the observer should 
see the source increasing in overall size and changing its shape from a full disk to a relatively 
narrow ring, at least while the expansion is relativistic. This is of importance for the possible 
gravitational microlensing of afterglows (Loeb & Perna 1997). If bolometric observations are 
obtained by piecing together band observations spanning many orders of magnitude, then most of 
the energy of the afterglow should be seen coming from a relatively narrow ring, at any time. 

The thickness of the zone that radiates most of the energy is important because it determines 
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the spread 5T in the arrival time of photons emitted at lab-frame t. Generally, 6T is comparable 
to T. In the radiative case electrons cool on a time-scale much shorter than the expansion time, 
only a very thin zone located behind the blast wave front releases significant energy and the 
thickness of the radiating fluid can be neglected. Therefore in this case the emitting fluid can 
be safely approximated as a surface. The eff'ect of the shell thickness is important only in the 
adiabatic case, and was taken in consideration by Waxman (1997b). We have employed here the 
equal arrival time surfaces using the kinematics of the blast wave. Radiation emitted by the fluid 
behind this shock is received at time T > Tsh, where Tg^ is the arrival time from the blast wave. 
Taking into account that the source size is increasing in time, it results that a flnite thickness leads 
to wider rings and lower values of the averages x and y given in Tables 1 and 2. Our estimations 
of the ring's width for an adiabatic remnant are larger by a factor up to 2.3 than calculated by 
Waxman (1997b). 

Additional complications arise if different regions on the equal-T surface are in different 
dynamic and/or radiative efficiency regimes. As the fireball decelerates, the cooling time-scale of 
electrons increases and they eventually become adiabatic. If there is a strong coupling between 
electrons and protons -|- magnetic field, the remnant and electrons evolve together from the 
radiative to the adiabatic regime. In the likely case that the coupling is weak (e.g. Meszaros , 
Rees, &: Wijers, 1997), the remnant evolution is adiabatic throughout, and only the electrons 
evolve from radiative to adiabatic. At times ^ 1 day, the most likely case is given by the third 
(or sixth) line of Tables 1 and 2. The real situation may be even more complex if the power-law 
electrons are in different radiative regimes, e.g. low energy electrons may be adiabatic {t'^y > t'^^p) 
while higher energy electrons may be radiative (t' < t' ). 
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TABLE 1. 

Intensity-averaged parameters on the equal arriving time surface 
and the width w of the ring seen by the observer for bolometric observations 



n [a, S] 


2(2n+l) 


^max 


^/(27oCT) 


y/ymax 


y/iiocT) 


Up/Uo 


w 


3.0 [0, 0]'' 


14 


0.78 


11 


0.87 


5.2 


39 


0.07 


1.5 [0, 1]'^ 


8 


0.82 


6.6 


0.76 


3.2 


4.0 


0.17 


1.5 [0,1]" 


8 


0.72 


5.8 


0.82 


3.4 


8.3 


0.11 


1.0 [2,0]'^ 


6 


0.74 


4.4 


0.77 


2.6 


6.7 


0.16 


0.5 [2, l]*^ 


4 


0.78 


3.1 


0.71 


1.8 


2.0 


0.23 


0.5 [2, 11" 


-J 


0. 19 


2.0 


0.76 


1.9 


G.8 


0.20 



^ radiative electrons " adiabatic electrons 



TABLE 2. 

Intervals for intensity-averaged parameters on the equal arriving time surface 
and width w of the ring seen by the observer for band observations 



n [a,d] 


^max 


x/(272cr) 


y/ymax 


y/iiocT) 


w 


3.0 [0,0]'^ 


0.88-0.72 


12-10 


0.77-0.89 


4.6-5.4 


0.16-0.06 


1.5 [0,1]'- 


0.87-0.78 


7.0-6.3 


0.70-0.79 


2.9-3.3 


0.25-0.13 


1.5 [0, 1]« 


0.89-0.77 


7.1-6.2 


0.67-0.79 


2.8-3.3 


0.30-0.14 


1.0 [2,0]'^ 


0.83-0.68 


5.0-4.1 


0.72-0.80 


2.4-2.7 


0.23-0.13 


0.5 [2, 1]'^ 


0.83-0.75 


3.3-3.0 


0.68-0.73 


1.7-1.8 


0.28-0.21 


0.5 [2.1]" 


0.78 0.5G 


3.1 2.3 


0.70 0.76 


1.7 1.9 


0.25 0.20 



^ radiative electrons " adiabatic electrons 
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Fig. 1. — Surfaces of equal arrival times, for a homogeneous external medium and a radiative 

(top) or adiabatic remnant dynamics (bottom). Each curve is a transverse section through the 
3-dimensional equal-T surface, highlighting the regions that radiate 50% (upper half of each 
curve) and 80% (lower half) of the bolometric flux. Projected on the plane perpendicular to the 
l.o.s. toward the center of explosion, these regions appear narrower than the projection of the entire 
radiating surface. The Cartesian coordinates are normalized to Tc^ecj which for the putative burst 
parameters in the text is ~ 4 x 10^^ cm, corresponding at a redshift z = 1 {Hq = 75kms^^Mpc^^, 
f2 = 1) to an angular scale 2.5 fias. 
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Fig. 2. — Distribution of peak of synchrotron spectrum on the equal arriving time surface (upper 
graphs) and bolometric luminosity on the same surface (lower graphs). Left panels are for a radiative 
remnant (n = 3), right panels for an adiabatic one (n = 1.5). The ring is shown by the steep rise 
in integrated bolometric luminosity in the lower graphs. Also shown in the lower graphs are the 
widths as observed in two fixed frequency bands (2 eV and 10 GHz). The lower right graph shows 
that if electrons are adiabatic (solid curve) , then the bright zone shrinks to an even narrower ring. 



